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Equilibrium and motion of a contact line are viewed as analogs of phase equilibrium and motion of an in- 
terphase boundary. This point of view makes evident the tendency to minimization of the length of the contact 
line at equilibrium. The concept of line tension is, however, of limited applicability, in view of a qualitatively 
different relaxation response of the contact line, compared to a two-dimensional curve. Both the analogy and 
qualitative distinction extend to a non-equilibrium situation arising due to coupling with reversible substrate 
modification. Under these conditions, the contact line may suffer a variety of chemo-capillary instabilities 
(fingering, traveling and oscillatory), similar to those of dissipative structures in nonlinear non-equilibrium sys- 
tems. The preference order of the various instabilities changes, however, significantly due to a different way the 
interfacial curvature is relaxed. 



I. INTRODUCTION 



Instabilities of a contact line are commonplace in various 
settings (for a recent review, see Refs. EE))- In many prac- 
tical applications, instabilities are detrimental when they, for 
example, impair smooth coating or spotless dewetting. Insta- 
bilities may lead, however, to formation of fine patterns or 
enhanced spreading in microfluidic applications, and cause 
such fascinating phenomena as spontaneous motion in non- 
equilibrium systems. Most thoroughly studied instabilities 
are of hydrodynamic origin, being caused, for example, by 
enhanced gravitational or Marangoni driving at a perturbed 
contact line, and are most commonly related to formation of 
capillary ridges [2|. Quantitative characterization of these in- 
stabilities is commonly based on lubrication approximation, 
and is strongly affected by the notorious contact line singular- 
ity resolved on a microscopic scale. 

We will concentrate here on a different kind of instabil- 
ity, caused by substrate modification rather than hydrodynam- 
ics, and therefore apt to occur at slow velocities and small 
scales. We start in Sect. [Iljwe discussing an analogy between 
equilibrium of a fluid body bounded by the contact line and 
thermodynamic phase equilibrium in two dimensions, which 
is made evident by variational representation of lubrication 
equations accounting for the action of surface tension and dis- 
joining pressure. We will argue, however, that this analogy 
does not justify extending the notion of two-dimensional line 
tension to the contact line, in view of a qualitatively different 
quasi-elastic response of the contact line dependent on three- 
dimensional effects hidden in the lubrication description. 

In the main part of the paper, we will develop the thermody- 
namic analogy in a different direction by considering chemo- 
capillary instabilities in a non-equilibrium setting arising due 
to coupling with reversible substrate modification. As a repre- 
sentative example, we will consider a set-up used in observa- 
tions of spontaneous motion of droplets induced by surfactant 
adsorption [3|. This system has been described by a model 
combining hydrodynamics in lubrication approximation with 
a linear reaction-diffusion equation for an adsorbed species, 
and studied both numerically (4) and analytically [5|. We re- 
consider it here as an analog of nonlinear models responsible 
for the formation of non-equilibrium structures 1 6 1 . There are 



remarkable similarities here, leading to the same variety of 
instabilities, which include, in addition to traveling instabil- 
ity, fingering instability and oscillations commonly found in 
other physical settings. The character of transitions is, how- 
ever, strongly influenced by the specific wavenumber depen- 
dence of the contact line relaxation response, different from 
the elastic response of an interphase boundary in two dimen- 
sions. 

We will derive equations of motion of a curved contact line 



in Sect. Ill by combining the effective equation of motion of 



a rectilinear contact line derived through resolving the contact 
line singularity [7 1 with the curvature response derived here in 
a local approximation and fitting the established theory 0][8]. 
After the diffusive surfactant model is introduced in Sect. |IV] 
the various instabilities of a contact line delimiting a static or 
moving semi-infinite fluid layer or bounding a circular droplet 
are analyzed in the following Sections. 



II. EQUILIBRIUM 

The evolution equation of the layer thickness h in lubrica- 
tion approximation can be written in a variational form |9] 



d t h = -V ■ n{h)V 



SW 



(1) 



The driving potential W is of the Cahn-Hilliard type, appro- 
priate for the case when the "order parameter" is conserved: 



W 
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The two terms in the integrand are the interfacial energy 
with the surface tension 7 and the net adhesion energy 
V is the two-dimensional gradient operator. Other extrinsic 
interaction terms, such as gravity, having the same algebraic 
form, can be added here, but will be further assumed to be 
negligible, as they operate on far larger scales. The mobility 
coefficient n(h) depends on h in a strongly nonlinear fash- 
ion when it is of hydrodynamic origin, e.g. n(h) = ^rp 1 /) 3 
for the Stokes flow with the dynamic viscosity r\ and no slip 
on the solid substrate. The validity of the lubrication hydro- 
dynamic approximation used in derivation of Eq. ([T]), as well 



2 



as of the concept of shape-independent adhesion energy, is 
formally restricted to films with a large aspect ratio, imply- 
ing small contact angles. It is universally used nevertheless 
for the analysis of both equilibrium and non-equilibrium films 
leading to qualiatively correct results even outside its formal 
applicability limits. 

It is clear from the variational form of Eq. ([T]i that the func- 
tional W is minimized in a stationary configuration, which is 
defined by the Euler-Lagrange equation 



7 v 2 /i-n(/i) = o, 



(3) 



where IT(/i) = $'(/i) is the disjoining pressure. One can con- 
sider two basic configurations. A macroscopic system (or a 
part thereof) can be viewed as an infinite layer asymptotically 
flat at spatial infinity. Practically, the layer may be shaped 
at distances far exceeding the range of adhesion forces by 
gravity or other macroscopic forces. The enormous macro- 
scopic to microscopic scale ratio, which may reach 7-8 or- 
ders of magnitude, leaves ample space for an asymptotically 
flat region at intermediate distances. The asymptotic slope is 
identified then with the equilibrium contact angle. 

Consider, for example, a commonly used expression for 
disjoining pressure of the type 



IL(h) 



A 



1 



(4) 



characterizing a partially wetting liquid with the net van der 
Waals interaction parameter A forming a precursor film of the 
thickness hf. The equilibrium contact angle 8 e is computed 
by solving the one-dimensional version of Eq. ([3]), 

jh"(x) -11(h) =0, h(-oa)=h f , ti(oo) = 6 e , (5) 

where x is the coordinate normal to the contact line. Using 
the phase plane representation h'(x) = p(h) and integrating 
yields 



[2*(M/7] 
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(6) 



A representative example of a disjoining pressure of a 
more complex form iflOl combines wetting van der Waals and 
dewetting polar interactions: 



n(/i) = -A + jBe -Vf 

h 6 



(7) 



In a dimensionless form with h scaled by the attenuation scale 
I and IT scaled by A/£ 3 , this function depends on a single 
parameter \ = B£ 3 /A: 



U(h) 



(8) 



The condition f(hf) = 0, which has positive roots at \ > 
Xf = (e/3) 3 « 0.7439, determines the thickness hf of a thin 
film in equilibrium with a flat macroscopic liquid layer. The 
equilibrium contact angle computed in the same way as ([6]) is 
positive at x > Xo = e 2 /8 0.9236. 

An alternative configuration is a microdroplet of a lim- 
ited volume (or several microdroplets). The volume con- 
straint can be treated by extracting from Eq. |2]i the term 



pV = p J h(x) d 2 x, where the Lagrange multiplier p plays 
a role of chemical potential. For the function Q, the simplest 
solution of the respective Euler-Lagrange equation, 



7 V 2 /i - 11(h) + p = 0, 



(9) 



is a parabolic cap of a radius R on the top of a precursor layer 
with the thickness somewhat exceeding hf. The Lagrange 
multiplier p can be identified then with the chemical potential 
or the vapor phase in equilibrium with the droplet of curvature 
radius R/6 e , i.e. with the vapor concentration just below the 
dew point when R^> hf9 e . For the function (|8jl, the solutions 
at at x < Xo an d a suitable value of p > can take the form of 
a straight-line front separating domains with alternative values 
of the stationary film thickness or of a circular "pancake" or 
"hole". 

Both macroscopic and microscopic configurations must be 
stable when the length of the contact line is minimal. Equa- 
tion |5]) or (j9j) with a suitable function H(h) can be recognized 
as equations of phase equilibrium with the layer thickness h 
playing the role of the order parameter. Respectively, Eq. ([TJ, 
is recognized as the Cahn-Hilliard equation governing relax- 
ation to the equilibrium state. The only unconventional fea- 
tures are a strongly nonlinear dependence of the mobility on 
the order parameter and a possibility for the latter to be infi- 
nite in a system of infinite extent. This analogy allows us to 
consider the question of contact line stability in the general 
framework of stability of phase equilibria. 

The question of stability of different configurations of an 
equilibrium contact line has been discussed in connection with 
the notion of line tension fffl . Unlike a line in two dimen- 
sions, deformations of a contact line inadvertently cause de- 
formations of the surface of the enclosed fluid volume, and 
therefore instability may not arise even when the line tension 
is formally negative [12]. Moreover, line tension of a contact 
line, unlike surface tension in three or line tension in two di- 
mensions, is not a fundamental property. If it is defined as a 
constituent part of the overall energy proportional to the line 
length, it remains dependent on the shape of the interface in 
the vicinity of the contact line, and therefore cannot be sepa- 
rated from surface tension and wetting properties. As noted in 
Ref. lfT2ll . stability of an equilibrium contact line is fully deter- 
mined by the applicable dependence of the disjoining pressure 
on the local layer thickness h, which determines both the equi- 
librium contact angle and the interface shape in the vicinity of 
the contact line affecting the apparent line tension. Mechkov 
et al |[T2l proved stability of a circular contact line and stabil- 
ity of a capillary ridge on short wavelengths. In the lubrication 
approximation, this result is evident: both the transitional area 
near the contact line and a front separating domains with dif- 
ferent film thickness carry positive energy, and their length 
should be minimized by the globally stable configuration. 



III. EQUATION OF MOTION FOR THE CONTACT LINE 

Instabilities of a contact line or a front between alternative 
equilibrium thickness levels may arise in a one-component 
fluid only in a non-equilibrium setting. The key ingredient 



3 



of the stability analysis in different situations is the equation 
of motion for the contact line. A naive approach based on the 
local curvature of the contact line and line tension that would 
be sufficient if the problem was truly two-dimensional is not 
applicable here, since the three-dimensional configuration of 
the interface in the vicinity of the contact line plays a cru- 
cial role. Strictly speaking, the problem is nonlocal, since any 
distortion of the contact line causes pressure changes within 
the liquid layer or droplet that affect all other locations along 
the contact line. The problem becomes tractable only when 
pressure fluctuations are neglected and the analysis is based 
on a local equation of contact line motion. Such an equation 
can be extracted from the analysis of slow displacement or 
spreading of droplets based on the solvability condition for 
perturbed stationary configurations J5J0Q2). The theory is 
based on matching of perturbations induced by displacement 
of the contact line in the macroscopic domain and in the con- 
tact line vicinity where the dynamic contact line singularity 
is resolved through either slip or the presence of a precursor 
film. For a droplet spreading due to a difference between the 
apparent macroscopic angle 9 and the equilibrium contact an- 
gle 6 e , the spreading velocity U is Q 



u = l -Qu {e z 



5 ). 



1 



3ln(aR/iy 



(10) 



where Uq — 7/77 is the characteristic viscous velocity and the 
argument of the logarithm contains the ratio of a macroscopic 
length R (droplet radius) to a microscopic length I (thickness 
of the precursor layer or slip length); the numerical factor a 
is dependent on the functional form of the disjoining pres- 
sure and/or other microscopic and macroscopic inputs. The 
linearized form of Eq. ( fT0] > applicable when the difference 
■d = 9 — 9 P is small is 



U 



E7b0. 



(11) 



If this expression is applied to compute the displacement ve- 
locity of a droplet under the action of a variable equilibrium 
contact angle, the result coincides precisely with that obtained 
earlier in a somewhat different way [5 |. We will further adopt 



Eq. ( 10 1 and its linearization as the local equation of motion 
of the contact line also for non-circular geometry, neglecting 
nonlocal effects which are likely to manifest themselves only 
as weak corrections of the numerical geometric factor in the 
argument of the logarithm. 

Perturbations of the apparent contact angle i? can be com- 
puted in the following way. Let the contact line be shifted 
from an unperturbed position T(y) along the outer normal by 
an increment £(y), where y is a spanwise coordinate. The 
displacement is supposed to be small on the macroscopic 
scale, although it may be large compared to the microscopic 
scale £; the essential requirement is that the spanwise deriva- 
tive £'(y) be small, so that the corrugation wavelength be 
much larger than t. Projecting the perturbation upon T im- 
plies the perturbation of the local film thickness h(T,y) — 
e £(y). The global perturbation is obtained by solving the 
two-dimensional Laplace equation V 2 /i = in the unper- 
turbed geometry, and d is subsequently computed as the nor- 
mal derivative n ■ V/i = at the contact line. 



Consider first relaxation of a perturbed rectilinear contact 
line bounding the liquid layer at x < with a constant incline 
9 e . The contact line perturbation is presented as a Fourier 
integral £(y) = J (k cos ky dk. For small deviations of the 
contact angle, the perturbation of the film thickness at the axis 
x = (i.e. at the nominal position of the undisturbed con- 
tact line) induced by a Fourier component with k <C t~ x is 
h(0, y) — e Cfc cos ky. The perturbation of the film thickness 
h(x, y) in the bulk of the film obtained by solving the Laplace 
equation with this boundary condition is 



h(x,y) = 9 e ( k e Wx cos ky. 
This yields the perturbation of the contact angle 

i? = -h x (0, y) = -\k\9eCk cos ky. 



(12) 



(13) 



Using this in ( 1 1 1 yields the Fourier component of the contact 
line displacement velocity 



Cit = -Qu e 3 e \k\a 



(14) 



We see that the restoring force is proportional to \k\ rather 
than k 2 as it would be for a truly two-dimensional elastic line. 
This is in agreement with the qualitative argument ll8l [T5l [T6ll 
stipulating that the capillary energy (ostensibly proportional 
to k ) should be integrated over the penetration length of the 
order |fc| _1 , thus leading to a reduced power of the wavenum- 
ber dependence. This makes the response of the contact line 
"superdiffusive", e.g. the curvature radius r of a small bulge 
grows with time as r a t 2 ^ 3 rather than r oc t 1 ! 2 . The re- 
laxation rate proportional to |fc| is also supported by full nu- 
merical computations [14|, which have shown that it breaks 
down only near the Landau-Levich entrainment limit. In this 
limit, the approximation breaks down, as the contact line per- 
turbation computed as above becomes comparable with the 
unperturbed contact angle. 

Another example to be studied in detail below is a circular 
droplet of the radius R. The unperturbed profile is a parabolic 
cap 



ho(r) 



9 e R 
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(15) 



where r is the radial coordinate r. The perturbation modes 
dependent on the polar angle <f> are £„ = £ n cos ncj) with an in- 
teger n. The solution of the Laplace equation for the perturbed 
interface is 



h{r,<j>) = 9 e ( n (r/R)' 



cos no 



(16) 



Taking also into account the contribution of h' (r) yields the 
perturbation of the apparent contact angle at the contact line 



$ = 9 e -ti {R + l; n )-h r {R,4>) = - 



n - 1 
R 



9 e Q n cos no 



(17) 

In a special case of a symmetric perturbation (n = 0), the 
perturbation of 9 is related to the perturbation of R by the 
conservation condition of the droplet volume V = ^9 e R 3 , 



which yields ■& = —( n 9 e /R. The expression ( 17 1 can be ex- 
tended also to this case if n — 1 is replaced by its absolute 
value In — 1|. 
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IV. DYNAMICS OF A DIFFUSIVE SURFACTANT 

A disturbance displacing the contact line may be caused 
by chemical inhomogeneities, which can be accounted for by 
assuming a linear dependence of the equilibrium contact angle 
on the local concentration of an adsorbed chemical: 



9o[l-6(c-c )], 



(18) 



where 6*o is the reference contact angle at c = Co and b is a 
proportionality constant. Modification of wetting properties 
by adsorption or chemical reactions has been used in a num- 
ber of experiments to induce spontaneous motion of droplets 
on solid substrate. In the following, we will consider a re- 
versible setup allowing for restoration of substrate properties 
0. Variation of the contact angle is caused in this system by 
deposition of a surfactant from a bulk phase (assumed to have 
negligible viscosity) and its dissolution underneath a droplet 
or liquid layer. The theoretical model describing stationary 
droplet motion in this setup [5 1 is based on a surface diffusion 
equation with a constant adsorption rate and linear desorption 
kinetics. This leads to the surfactant diffusion equation writ- 
ten in the dimensionless form 



d = V 2 c- c + H(x). 



(19) 



The surfactant coverage c is scaled here by the coverage in 
equilibrium with the constant surfactant concentration in the 
continuous phase, time by the inverse desorption rate constant 
kd, and length, by y/D/kd, where D is the surfactant diffusiv- 
ity on the substrate; H (x) is the Heaviside step function equal 
to 1 outside and inside the footprint of the droplet or liquid 
layer. For simplicity, it is assumed that the desorption rate is 
constant all over the substrate and the bulk concentration in 
the liquid layer is negligible. 

We will apply now the same model in conjunction with 



Eq. ( 10 1 to explore stability of a contact line. In this formu- 
lation, the problem is almost identical to that of stability of 
non-equilibrium structures in the FitzHugh-Nagumo or simi- 
lar models (6|. First, we have to find the stationary solution of 
Eq. 



19 1 for the unperturbed configuration. For a semi-infinte 
layer at x < 0, the solution dependent on a single coordinate 

x is 



at x < 0, 
at x > 0. 



(20) 



The reference concentration Co can be identified with the sta- 
tionary concentration at the contact line c s (0) = 1/2. 

For a circular droplet, the stationary solution c = c s (r) of 



Eq. (19 1 is 



= J RKi(R)I (r) at 
< 1 - Rh(R)K (r)} at 



r <R 
r>R, 



(21) 



where I n , K n are modified Bessel functions. The values of 
c s (R) given by both lines of this formula are the same, in 
view of the identity Ki(R)I (R) + K (R)h(R) = 1/R. 
This common value should be identified with the reference 
concentration cq. 



The concentration perturbation due to a contact line dis- 
placement can be computed in a simple way by observing that 
any shift of the contact line position by an increment £ <C 1 is 



equivalent to switching on or off the source term in Eq. ( 19 1 
in a narrow region near the contact line. This contributes to 
the equation for the perturbation c = c — c s a source localized 
at the unperturbed contact line position and proportional to its 
displacement: 



d t c = V 2 c-c~Z5(T). 



(22) 



The dynamic equation for the contact line displacement is 
obtained by linearizing Eq. ( 10 1. Using Eq. (fl"8]l and expand- 



ing the surfactant concentration in the vicinity of the unper- 
turbed contact line position yields 6 e = 9q[1 — b(c + 
where j is the stationary surfactant flux into the liquid across 



the contact line. Combining this with Eq. (Hi yields the equa 



tion for the local contact line displacement 

£ = -xWOo-b(c + jQ]. 



(23) 



We have transformed here to the dimensionless form based 
on the chemical length and time scales. The remaining di- 
mensionless parameter x = QUoOg/y/Dkd is proportional to 
ratio of the characteristic capillary velocity Uq ~ j/i] to the 
characteristic chemical velocity \jDkd- This ratio is typically 
very large, but is effectively reduced being multiplied by the 
cube of the contact angle, which must be small when the lu- 
brication approximation is applicable. 

Although perturbations are habitually presumed to be ar- 
bitrarily small for the purpose of linear stability analysis, the 
method is actually also applicable to finite perturbations of the 
contact line position, provided they are small compared to the 
characteristic surfactant diffusion scale. 



V. INSTABILITY OF A RECTILINEAR CONTACT LINE 

For a rectilinear contact line displaced by a harmonic per- 
turbation, d in Eq. ( 23 1 is given by Eq. ( pj) and the station- 
ary flux following from Eq. (20i is j = c' s (0) = 1/2. Pre- 



senting the concentration perturbation in the spectral form 
c(x,y,t) = rp(x,t) cos ky and setting in the perturbation 



equations ( 22 1, ( 23 1 ip, Cfe 
lem 



e At leads to the eigenvalue prob- 



Kk = -xMCk-bty + Ck/2)], 
ipxx - q 2 ip = (kS(x), 



(24) 
(25) 



where q 2 = 1 + A + k 2 . The solution of Eq. (25 I, presuming 
Re q > 0, is 



V>(x) 



2q 



-q\x\ 



Using this in Eq. (|24| yields the dispersion relation 



A = -x 









H)] 







(26) 



(27) 
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FIG. 1 : Monotonic and oscillatory instability regions in the paramet- 
ric plane x, b for a rectilinear contact line. 



The stability condition is Re A < 0. As expected, A van- 
ishes at k = 0, which reflects the translational symmetry of 
the rectilinear line. Unlike a similar problem in Ref. O with 
fc replaced by fc 2 , the derivative dA/dfc at k = is always 
negative, so that the contact line is stable to long-scale pertur- 
bations. Stability is lost as A vanishes at 



'l + V5 



(3 + \/5) 3 / 2 



(28) 



A band of unstable modes around k = k c (which turns out to 
coincide with the golden ratio) widens at b > b c . In dimen- 
sional units, the wavelength of growing perturbations is of the 
same order of magnitude as the diffusional length y/D/kd- 

In addition, oscillatory instability is possible when \ is suf- 
ficiently large. It first arises in the long-scale low-frequency 
mode at x = Xc = 4/6, and can be detected by expanding 
Eq. (27 1 near this point. For x — Xc = c *C 1, the suitable 
scaling is 



A = eA,, 



ie 1 / 2 ^ 



k = ek. 



(29) 



The leading terms of the real and imaginary parts of the ex- 
pansion are, respectively, of the order 1 and 3/2. This yields 
two lowest order equations for the real and imaginary parts of 
A as functions of fc: 



fc bZ) 2 = 0, 

16 



6w(26-12A r -5w 2 ) =0. (30) 



The first of these gives the frequency as the function of the 
wavenumber: 




(31) 



The second equation gives the O(e) real part of the eigenvalue 



20 fc 
~9b' 



(32) 



indicating that instability indeed occurs at x > Xc an d dis- 
appears as fc increases. Higher orders of the expansion give 



0(e 2 ) corrections to the real, and 0(e 3 / 2 ) corrections to the 
imaginary parts of A. Within the instability region, wave 
modes with finite wavelength also become unstable and are 
likely to exhibit there the fastest growth rate. The monotonic 
and oscillatory instability regions in the parametric plane x, b 
are shown in Fig.[T] 

VI. INSTABILITY OF A CIRCULAR DROPLET 



For a circular droplet, $ in Eq. (23 i is given by Eq. ( 17 1 



and the eigenvalue problem for coupled spectral equations of 
the displacement and concentration perturbation analogous to 



Eqs. (24 1, (25 1 has the common form for any integer n > 0: 



AC, 



-x[|n-l|C„-6(V+Kn)], 

,2 



(33) 



^ + -V^V - I <? + ^- ) *l> = Cn5{r - R), (34) 



where q = 1 + A. The solution of Eq. ( 34 1 is 



-R( n K n (qR)I n (qr) 
-R( n I n (qR)K n {qr) 



r<R, 
r > R. 



(35) 



Using this, together with the expression for the stationary flux 
j = cJ s (R) = RKi(R)h(R) following from Eq. ((21}, in ((33 1 
yields the dispersion relation 



A 

X 



II 



R 



bRlhWKiiR) - I n {qR)K n {qR)] 



(36) 

As expected, A vanishes at n = 1, which corresponds to a 
shift without deformation. For n^l, the critical value b = b n 
for a given radius verifies Eq. ([36| with A = 0: 



1 



R 2 



II 



[I 1 (R)K 1 (R)-I n (qR)K n (qR)}. (37) 



For n = 0, the value bo given by this relation is negative, so 
that no monotonic instability can arise in the symmetric mode. 
The most dangerous instability mode is the one with the low- 
est critical value of b±(R) — min n b n (R). For small droplets, 
the dipole mode deforming the circular droplet into an ellipse, 
is the most dangerous one. With the growing radius, the most 
dangerous monotonic instability mode shifts to larger values 
of n, as shown in Fig. [2] The absolute lower limit of mono- 
tonic instability is 62 ~ 5.158 at R sw 2.337, which is below 
the respective limit b c » 6.6604 for the rectilinear contact line 
given by Eq. ( |28j ); the latter limit is approached asymptotically 
at R —> 00, as the envelope of the family of curves in Fig. [2] 
gradually rises. The wavenumber fc = n/R rises as well on 
the average, as seen in the inset in the same Figure, coming 
gradually closer to the value fc c 1.272 of Eq. (28 1. 

Besides the monotonic instability, dynamic instabilities de- 
pendent on the parameter x are possible. Those are traveling 
instability setting the droplet into motion and oscillatory or 
wave instability. The threshold of traveling instability is de- 
tected, according to the general algorithm f6|, by differenti- 
ating the dispersion relation ( 36 1 for n = 1 with respect to A 



6 




FIG. 2: Monotonic instability thresholds for a circular droplet of ra- 
dius R. The numbers indicate the value of n. Instability occurs above 
the lower envelope of the family of curves. Inset: the wavenumber 
k — n/R for the most dangerous mode; the dashed line indicates the 
respective value for the rectilinear contact line. 
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FIG. 3: The lower limit of the product b\ at the traveling instability 
limit as a function of radius R. Inset: The lower limit \t at which 
the traveling instability precedes the monotonic one. 



and solving the resulting equation at A 
into spontaneous motion at 



0. The droplet sets 



R[Iq(R)Ki(R) - h{R)K 2 {R)] > 2/b X . 



(38) 



The lowest value of the product b\ enabling traveling insta- 
bility bx ~ 11.46 is attained at R w 1.5866 (see Fig. [3}. The 
lower limit \ — Xt at which the traveling instability precedes 
the monotonic one for droplets with a certain radius can be ob- 
tained by using in Eq. (38 i the critical value b± for the most 



dangerous mode given by Eq. (37 1. The plot Xt{R) is shown 
in the inset of Fig. [3] 

Oscillatory instability can be detected numerically by solv- 
ing coupled equations for the real and imaginary parts of 



Eq. ( 36 1. Setting A = \uj we obtain at the instability threshold 

1 R 2 



n-l 



[h(R)Ki(R)--RoF n {R,u)], (39) 



- = -— ImF n (R,u>), 
X w 



(40) 



F n (R,oj) = I n (RVl + kJ) K n (RVl + iu) ■ (41) 



Using Eq. ( 39 1 to eliminate b, we are left with a single equa- 
tion d40} , which can be solved numerically to compute the 
frequency u> at the Hopf bifurcation point for a droplet of a 
given radius R. 

For n > 2, the locus of oscillatory instability branches off 



the locus of monotonic instability defined by Eq. ( 37 1 at the 
point of double zero eigenvalue where the frequency vanishes. 
Close to this line, u is small, and a Taylor series can be used. 



The leading term in the expansion of Eq. ( 39 1 is of the order 
0(w 2 ), while the leading term in the expansion of Eq. (40 1 
is linear in ui. Setting the latter to zero gives the value of 
X = Xn required to obtain the double zero as a function of 
the radius R; the same relation can be obtained most readily 



by differentiating Eq. (36i with respect to A. The resulting 



(42) 



condition is a generalization of Eq. 

R[I n ^(R)K n (R) - I n (R)K n+1 (R)] > 2/b Xr , 



The value of \n given by this relation consistently rises with 
n, and is always higher than the respective value for traveling 
instability obtained for n = 1. This proves that low-frequency 
oscillatory instability never occurs before traveling instability 
sets in. 

The remaining possibility is oscillatory instability with a 
finite frequency, which is also feasible (and, indeed, is most 
likely to occur) in a symmetric mode. Rather than solving 



Eqs. (39 1, (40 1 directly, it is more instructive again to com- 
pute the lower limit \ — \o at which this instability precedes 
the monotonic one for a droplet with a certain radius. For 
this purpose, uj is computed by solving numerically Eq. ( f39| ) 
with b = b±, and then \o is evaluated using Eq. (40 1. The 



computation shows that the frequency at the Hopf bifurcation 
decreases monotonically, going down to zero as the limit of 
a rectilinear contact line approaches, in agreement with the 
results of Sect. [V] The curve Xo(R) for the symmetric mode 
lies, however, slightly higher than the curve xt{R), so that 
traveling instability always occurs first. 



VII. INSTABILITY OF A MOVING CONTACT LINE 

The analysis can be extended to the case when the differ- 
ence between the equilibrium and apparent contact angles is 
not small and the unperturbed line is moving with a certain ve- 
locity U. The deviation d is now redefined as a perturbation of 
the apparent contact angle at the moving contact line 9, which, 
in accordance with Eq. ( 10 1, equals to (9 z e + 3U/QU ) 1/3 . In 

[ffl] $ is given for a rec- 
13 i with 8 e replaced by 9. The 



the same approximation as in Sect, 
tilinear contact line by Eq. 



equation of the local contact line displacement ( 23 1 has to be 
written now for the displacement £ relative to the unperturbed 
moving line, and is modified to 



i = -(x + 3u) #/e + b X (c+jt), 



(43) 



where u — U / \fkdD is the dimensionless velocity. 

The flux through the contact line defined by solving 
Eq. ( 19 1 in the comoving frame is j — c' s (0) = 1/V4 + u 2 , 
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FIG. 4: Loci of monotonic (dashed lines) and oscillatory (solid line) 
instability in the parametric plane x, b for the rectilinear contact line 
moving with the speed it = ±0.1 (as indicated by numbers at the 
curves). The stable domain is in below both respective curves. 



and the spectral equations ( 24 1, ( 25 i are modified to 

Ck 



Kk = -(x + 3u)|fc|Cfc + &x(^- 

ipxx - ui/j x - q 2 ip = ( k S(x) 
The solution of the last equation is 

0= 



V4" 



(44) 
(45) 



ip(x) = - 



Vm 2 



exp 



(46) 

where the positive and negative signs apply, respectively, at 
x > and x < 0. Using this in Eq. ( |44| yields the dispersion 
relation 



A= -{ X + iu)\k\+b X 



1 



1 



(47) 



Notably, the coefficient at |fc| vanishes at u = — x/3, which 
corresponds to the Landau-Levich entrainment limit; the the- 
ory becomes inapplicable close to this point fl4l . Besides 
the coefficient at |fe|, the dispersion relation is insensitive to 
reversing the direction of motion, but the reference concentra- 
tion Co defined as the surfactant concentration on the unper- 
turbed contact line, is different in the two cases. 
Monotonic instability is observed at 









u 2 \ 



















(3 + 75) 3/2 (l + U 2 /4) (3u + X ) 



(48) 
(49) 



xvi+71 

Besides the velocity dependence, a qualitatively significant 



change, compared to Eq. ( 28 1, is the dependence of the mono- 
tonic instability threshold on the parameter x which appears 
at any non-zero velocity; the marginal wavelength k c remains, 
however, independent of x- At x 0, the velocity depen- 
dence becomes singular: a retreating contact line u < is 



always unstable under these conditions, while for an advanc- 
ing line (u > 0) the threshold rises sharply. 

Oscillatory instability emerges, as at u = 0, in a long- 
wave low-frequency mode, and can be detected by expanding 
Eq. (47 1 using the scaling (29 1 as before. The 0(e 1 / 2 ) term 



gives the parametric relation at the instability threshold 

b x = 4(1 + w 2 /4) 3/2 . (50) 
The dispersion relation uj(k) following from the O(e) term is 



LO 



k (u 2 + 4) 



(u 2 + 4) 3/2 



6 u 



(51) 



For u > 0, this relation is qualitatively similar to Eq. ( f3T) ; 
the function \ r (k) defined by the 0(e 3 / 2 ) term is likewise 
similar to Eq. ( [32| . The overall bifurcation diagram changes, 
however significantly, compared to that in Fig. [T] The loci of 
monotonic and oscillatory bifurcation fail to intersect already 
at u > 0.2012; at higher velocities, only oscillatory instability 
is relevant. At u < 0, on the contrary, monotonic instability 
prevails at small x ( see Fig-Q. 

VIII. DISCUSSION 

We have seen in Sect.|II]a remarkable similarity between the 
contact line equilibrium (described in lubrication approxima- 
tion) and equilibrium of an interphase boundary in two dimen- 
sions. Likewise, chemo-capillary instabilities of the contact 
line driven out of equilibrium by coupling with surfactant ad- 
sorption (or other non-equilibrium process affecting wetting 
properties of the substrate) are remarkably similar to instabil- 
ities of two-dimensional non-equlibrium structures. A sub- 
stantial difference stems, however, from relaxation response 
of a curved contact line being proportional to the absolute 
value rather than square of the wavenumber k. This may be 
attributed to a "two-and-a-half-dimensional" character of the 
lubrication approximation. In view of this specific response, 
the concept of line tension lifted from the two-dimensional 
world is not duly applicable to the contact line. This is mani- 
fested, in particular, in the paradox of stability of the contact 
line at an apparently negative line tension lfTTl[T2l . 

Out of equilibrium, the modified wavenumber dependence 
brings about more subtle changes, compared to instabilities 
of dissipative structures in two dimensions HQ. The onset of 
instability of a rectilinear contact line shifts to a finite wave- 
length, and traveling instability prevails over the oscillatory 
one for circular droplets. We have seen in Sects. fV} | VI| that 
monotonic instability occurs on the wavelength of the same 
order of magnitude as the surfactant diffusion length. It should 
cause fingering in the case of macroscopic liquid volumes of 
a size far exceeding this length, like it is commonly observed 
in spreading of surfactant-laden liquid layers ifTTl . It can also 
explain blebbing and fission of moving droplets observed in 
the experiment [18|. Instabilities with lower n may result in 
splitting of microdroplets with a radius of the same order of 
magnitude as the diffusion length. Dynamic instabilities are 
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enhanced by the increased ratio of the characteristic hydrody- of instabilities predicted by theory may be masked in actual 
namic and reaction-diffusion velocities expressed by the pa- experiments by surface roughness arresting contact line dis- 
rameter \- One should be aware, however, that the variety placement. 
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